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C/^ Abstract. 

, — ^ We introduce a systematic framework to calculate the bipartite entanglement 

^"*, entropy of a compact spatial subsystem in a one-dimensional quantum gas which 

can be mapped into a noninteracting fcrmion system. We show that when working 
^_^ with a finite number of particles A'^, the Renyi entanglement entropies grow as 

InAf, with a prefactor that is given by the central charge. We apply this novel 

technique to the ground state and to excited states of periodic systems. We 

fl\ also consider systems with boundaries. We derive universal formulas for the 

e leading behavior and for subleading corrections to the scaling. The universality 

of the results allows us to make predictions for the finite-size scaling forms of the 
, 1^ corrections to the scaling. 






I 

o 
o 



(N 
> 

in 

00 

en 

o 



X 



The entanglement entropy of ID systems in continuous and homogenous space 2 

1. Introduction 

Entanglement is a fundamental phenomenon of quantum mechanics. Much theoretical 
work has focused on the entanglement properties of quantum many-body systems, 
showing their importance to characterize the many-body dynamics [T]. In particular, 
lots of studies have been devoted to quantify the nontrivial connections between 
different parts of an extended quantum system, by computing von Neumann or Renyi 
entanglement entropies of the reduced density matrix pA of a subsystem A. Renyi 
entanglement entropies are defined as 

S^ = -^\nTTpl. (1) 

1 — a 

For a — > 1 this definition gives the most commonly used von Neumann entropy 
5*1 = — trp^lnp^, while for a — > oo is the logarithm of the largest eigenvalue of 
Pa also known as single copy entanglement [2]- 

One of the most remarkable results is the universal behavior displayed by the 
entanglement entropy at ID conformal quantum critical points (i.e. with dynamical 
critical exponent z — 1), determined by the central charge [3] of the underlying 
conformal field theory (CFT) 2l[5l|6|[7]. For a partition of an infinite ID system 
into a finite piece A of length i. and the remainder, the Renyi entanglement entropies 
for (. much larger than the short-distance cutoff a are 

^a = |fl+-Vn-+Ca, (2) 

6 \ aj a 

where c is the central charge and Ca a non-universal constant. When ^ is a finite 
interval of length ^ in a finite periodic system of length L, CFT predicts the universal 
asymptotic scaling [B] 

5'a = ^ (l + -)ln(— sin^) -hca, (3) 

o\ a J \Tra L J 

where, remarkably Ca is the same non- universal constant in Eq. (pi). For future 
reference, it is also important to mention the result in a finite system of length L with 
some boundary conditions at its ends and for an interval of length H starting from one 
of the two boundaries [6j IH |9] : 

C /, 1 \ , /2L . -kIx Cn 



^„ = -^l + -^Jln(-sin-J + ^+ln,, (4) 

where again Cq is the same non- universal constant as above and Ing is the universal 
boundary entropy of Affleck and Ludwig [lU] . All the Renyi entropies Sa are proper 
and equivalent measure of entanglement in a pure state [T] , but the knowledge of Sa 
for different a characterizes the full spectrum of non-zero eigenvalues of pA (see e.g. 
[llj ) providing significantly more information on the entanglement than the solely 
knowledge of the von Neumann entropy. 

The CFT results reviewed above have been confirmed in many spin chains and 
in ID itinerant systems on the lattice (too many to be mentioned here, we remand 
the interested reader to the comprehensive reviews on the subject ^). These studies 
have allowed a deeper understanding of the convergence and precision [12] of ID 
simulation algorithms based on the so-called matrix product states [13]. However, 
analogous results must also be valid for systems in continuous space, and therefore 
directly derivable in continuous models. Apart from the interest to describe trapped 
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ID gases experimentally realized with cold atoms, the entanglement of continuous 
models is also instrumental to develop ID tensor network algorithms for gases, as the 
one proposed in [Hj. Despite of this fundamental interest, almost no effort (with the 
exception of Refs. [THIIII] and the orbital partitioning in quantum Hall states [T7]) has 
been devoted to the spatial entanglement of gas models (that is distinguished from 
the particle partitioning [18]). 

In a previous short communication il9] we introduced a systematic framework to 
tackle free fermion gases in any external conditions for an arbitrarily large number of 
particles. The most general result of this investigation was that, when dealing with 
a finite number of particles A'^, the ID entanglement entropy grows like In TV, with a 
prefactor that again is given by the central charge. In this formulation N acts as an 
UV cutoff, representing a concrete alternative to the lattice. In this manuscript, we 
detail the calculations in Ref. [T^ for homogenous ID gases and we report a series 
of new results about the leading and subleading corrections to their scaling behavior. 
The degree of universality of these results allows to us to make novel predictions 
for spin chains on some universal functions describing the corrections to the scaling. 
The determination of these functions was left as open problem by previous lattice 
investigations pOll^l^ . 

1.1. The model and its eguivalence with others 

We consider a system of free spinless fermions in the continuum interval [0, L]. We 
work with a finite number of particles N. Therefore all the quantities and in particular 
entanglement entropies are finite since N acts as a cutoff. Appropriate boundary 
conditions (BC) are imposed in order to have a discrete energy spectrum. 

Apart the per se interest, spinless free fermions are also equivalent to other models 
of direct physical application. The ID Bose gas with short-ranged repulsive interaction 
(i.e. the Lieb-Liniger model [23]) with Hamiltonian 

^^ = -Ea^+2C E '5(x,-xO, (5) 

j = l 3 l<j<l<N 

in the limit of strong interaction C — > cx) (i.e. impenetrable bosons, also known as 
Tonk- Girardeau gas) is exactly mapped to spinless fermions [24 and the entanglement 
entropy of a single interval in the two models do coincide, because the boson in an 
interval are functions only of the fermions in the same interval (this is not true anymore 
in the case of more disjoint intervals because of the presence of a bosonization string, 
analogously to spin-chain models [25i i26i |27j). The properties of the Lieb-Liniger 
model are described solely by the dimensionless parameter 7 = CL/N 23J, thus the 
Tonks-Girardeau limit describes the dilute model (i.e. N/L ^1) for any value of C. 
Another important model mappable to free fermions is the spin-1/2 XX chain 
defined by the Hamiltonian 

^ 1 
where a^'^'^ are the Pauli matrices at site /. The Jordan- Wigner transformation 

\m.<l J 
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maps this model to the quadratic Hamihonian of spinless fermions 

^ = - E ^^'+1 + 4iCi + 2/i {c\ci - J) . (8) 



2, 

Here h represents the chemical potential for the spinless fermions c/, which satisfy 
canonical anti-commutation relations {q, cJ„} = (5;,m- The Hamiltonian m\ is diagonal 
in momentum space and for |/i| < 1 the ground-state is a Fermi sea with filling 

^^^ arccos|fe| _ (9) 

TT 

Only for |/i| < 1 we are dealing with a gapless theory. 

The continuum limit of the Hamiltonian ([8]) is then the system of free-fermions 
we are considering in this paper and so all the universal properties that do not depend 
on lattice regularization can be obtained from the continuum model. At this point, 
it is worth discussing how to obtain the continuum limit in some details. The lattice 
model is formed by Liat sites separated by the lattice spacing a (usually set to 1 in 
all lattice studies). N particles populate the chain with filling v ~ iV/Liat and we are 
interested in the entanglement entropy of £iat sites. The continuum limit is a system 
of A'^ free fermions in a box of length L and is obtained by sending a — > 0, j^ ^- 
keeping fixed i/^iat equal to iN/L in the continuum, where L = aLiat and (. = a£iat. 
This allows us to use the CFT results to predict a priori some of the results we are 
going to derive. In Ref. |2H], Eq. (pi) has been derived for the XX chain, obtaining 

^^^-g(l + ^)ln(^2sin7r^)+i?„, (10) 



1 \ f"^ dt 

^ 



(11) 



where also the non-universal constant is determined as 

1 / 1 _ 1 \ 1 _ e-2 

a^2 \asmht/a sinhtj sinht 6 

Combining this exact result with the CFT prediction in Eq. ([3]), we get the asymptotic 
scaling behavior of the entanglement entropies in finite XX chain 

Ea. ■ (12) 



S^^ = l(l + -]\n 



sm I 1 1 sm ixv 

na \ Liat 



Taking now the continuum limit a — > 0, i^ — > as explained above, we arrive at the 
prediction 






1 / 1 



6 V a 



If i\ . /... . e 



•^^ « 1 + - ln(l/a) -- 1+- In 2iVsin7r- +£:„, (13) 



6 V a J \ L 



where the subtraction of the term proportional to In a comes from the normalization 
of the reduced density matrix. 

Finally we quote the ID Bose-Hubbard model described by the Hamiltonian 

^BH = ^ 5](6.+i - b,)\b,+i -b,) + ^Yl '^'("^ - 1)' (14) 

■i i 

where bi are bosonic operators and n^ = 6^6^ is the particle density operator. The 
hard-core limit C/ — >■ 00 of the Bose-Hubbard model implies that the particle number 
Hi per site is restricted to the values nt = 0,1, and so in this limit can be exactly 
mapped into a lattice model of spinless fermions. Clearly the continuum limit of the 
Bose-Hubbard model is nothing but the Lieb-Liniger gas introduced above. 
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2. The method 

We consider a system of N non-interacting spinless fermions with discrete one-particle 
energy spectrum, such as a finite system or one confined by a proper external potential. 
The many body wave functions ^^(xi, ...,XAr) can be written in terms of the one- 
particle eigenstates as a Slater determinant 

*(xi,...,XAr) = — =det[(/)fc(a;„)], (15) 

where the normalized wave functions 0fc(a;) represent the occupied single-particle 
energy levels. The ground state is obtained by filling the N levels with lowest energies. 
Thus, the ground-state two-point correlator is 

N 

C(x, y) ^ {c\x)c{y)) = ^ (t>l{x)My) , (16) 

fc=i 
where c{x) is the fermionic annihilation operator and the one-particle eigenfunctions 
4>kix) are intended to be ordered according to their energies. The reduced density 
matrix of a subsystem A extending from xi to X2 can be written as 

PAOcexpf-/ dyidy2c\yi)'H{yi,y2)c{y2)\ , (17) 

J xi 

where % = ln[(l — C) /C] and the normalization constant is fixed requiring Trp^i = 1. 
This equation can be straightforwardly seen as the continuum limit of the formula 
for lattice free fermions [291 ISO], but can also be obtained following the standard 
derivation in Ref. [29] in path integral formalism. In the passage from lattice to 



continuum, the normalization factor in (17 1 depends explicitly on the lattice spacing 



a and it is responsible for the subtraction of the term proportional to In a in Eq. ( 13 ) 



We want to compute the bipartite Renyi entanglement entropies defined as in Eq. 
M of the space interval A in this fermion gas. For this purpose, we introduce the 
Fredholm determinant 

Da{X) = dct[X6A{x,y) - CA{x,y)] , (20) 

where Ca{x, y) is the restriction of C(x, y) to the part at hand from xi to X2 that can 
be written as CA{x,y) = 9{x — Xi)6{x2 — x)C(x,y)9{y — Xi)9{x2 — y) (or in matrix 
form Ca = PaCPa, where Pa is the projector on the interval A). The same definition 
holds for 6A{x,y) = PaS{x — y)PA- Following the ideas for the lattice model [5H], 
Da{X) can be introduced in such a way that it is a polynomial in A having as zeros 



the N eigenvalues of Ca- The Gaussian form of pa in Eq. (17) allows us to exploit 
the relation between the eigenvalues of pA and Ca to write 

^^^^^f|X^.nnD^^ (21) 

1 — a J Zin d\ 

X A Fredholm determinant is the extension of the standard determinant to continuous matrices. 
Its simplest operative definition is through the generalization to continuous kernels K{x, y) of the 
standard identity for determinants of a finite matrix M 

°° TrM^ °° TrK^ 

Indet [AI - M] = - ^ ^^ =^ Indet [\S(x ^ y) - K{x,y)] = - ^ -—^^ , (18) 

fc— 1 fc— 1 

where the traces are simply 

TrK" = dxidx2 . . ■dx„K{xi,X2)K(x2,X3) . . . K{x„-i,x„)K(xn,xi) . (19) 
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where the integration contour encircles the segment [0, 1], and 

ea(A) = -^ln[A" + (l-A)"] . (22) 

1 — a 



For a — > 1, ei(A) = —xhix — (1 — x)hi(l — x) and Eq. (21 1 reproduces the von 
Neumann definition. The integral representation (21) has been already derived and 
used in the context of discrete chain models [55] , thus involving the determinant of a 
standard matrix with the lattice sites as indices. 

The Fredholm determinant is turned into a standard one by introducing the NxN 
reduced overlap matrix A (also considered in Ref. |15| ) with elements 

/•X2 

A„™ = / dz(t>n{z)4>m(.z), n,m^l,...,N, (23) 

such that 

TrC^= / [lldy,)CA{yi,y2)-.-CA{yk-i,yk)CA{yk,yi) = 

rx2 k N N 

= / (n^y^ ^ rnM)^nAy2)-..Yl Kiy'^)^nAyi) = m 

= X]---H/ '^yiKAyi)(t^ndyi) ■ ■ ■ '^2/fc'/>^,(yfc)'^nfc-i(yfc) = TrA''", 

where from the first to the second line we use Eq. (16). Thus 

In Da{X) = - E ^ = - E ^a^ ^ l^det [AI - A] = ^ ln(A - a„) , (25) 

k—l k—1 jn—1 

where we use twice Eq. (flSl) and we denote with Um the eigenvalues of A. Inserting it 



into the integral (21), we obtain 



Sa{xi,X2) = 9 TT" E T^ ^ E ^aiara), (26) 

/ 27r2 ^ A - am ^ 

771 — 1 771 — 1 

as a consequence of the residue theorem. 

The matrix A is easily obtained for any non-interacting model from the one- 
particle wave functions, as the definition ( 23 ) shows. Calculating the entanglement 
entropies is then reduced to an A^ x TV eigenvalue problem that can be easily solved 
numerically and in some instances even analytically, as we are going to show. 

3. The ground-state of systems with periodic boundary conditions 

In a system of length L with periodic boundary conditions (BC), the normalized one- 
particle wave-functions are plane waves with integer wave numbers 

p27Tikx/ L 

Mx) = — 7^, keZ, (27) 

VL 

and energy Ek — 271^ k'^ / L? . For some physical problems, one has to impose anti- 
periodic BC for the fermion degrees of freedom (the appropriate BC can also depends 
on the parity of N) and so the momentum is quantized in terms of semi-integer wave 
numbers. However, as long as we are interested in the entanglement entropy of a 
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single interval, this does not change the final results because, as we shall see soon, the 
elements of the matrix A depends only on the difference between momenta that are 
always integer. Different results would be instead obtained for the entanglement of 
two disjoint intervals, similarly to what happens in CFT [31] |32] and lattice models 

E3[2ai2Zj- 

The element of the overlap matrix between two one-particle eigenstates with wave 
number fci and fc2 is 

L J XI 

7r(fci ~ fc2) 

The elements of the matrix A are not invariant under translation because of the 
explicit dependence on xi + X2 of the phase factor. However, the eigenvalues of A do 
not depend on this phase factor (but the eigenvectors do) and so also the entanglement 
entropies are translational invariant, as they must be. Indeed, in the determinant of 
AI — A, for each column we can bring out of the determinant the factor e^'^*'^i(^i+^2)/i 
and for each row e^'^^^'^^^^'^'^^ ^ ^ . Since fci and A;2 run on the same set of integers, the 
product of all these phases is 1, regardless of the values of the fcs. In the following we 
use this freedom to fix the phase factor to 1 and denote (. = X2 — xi. 

The ground state of a fermion gas with N particles is obtained by filling the N 
fc-modes with lowest energies. In the case of odd A'^, this amounts to fill symmetrically 
the A^ states with |fc| < {N — l)/2 (the zero mode is clearly included). For even N , 
there are two degenerate states obtained from the A^ — 1 ground state by filling the 
first available state either on the right or on the left. This small difference between 
odd and even terms does not play any role because the elements of the matrix A in 



Eq. (28) depend only on the differences between ki and ^2- Thus we can just start 



counting modes from the lowest k one-particle occupied state and the resulting matrix 



for a segment of length £ = X2 — a;i is from Eq. (28) 
svin:{n — m)(./ L 



, n,m = l,...,7V. (29) 

■K[n — m) 



By inserting the N eigenvalues of A into Eq. ( 26 1 , we obtain the entanglement entropy 



in a system of N particles. This is very easily done numerically as in Fig. [T] 

3.1. The leading behavior of the entanglement entropies 

We can also compute analytically the large N behavior of the entanglement entropies. 
Indeed we can write In D^ = In det G, where G = AI — A is an iV x A^ Toeplitz matrix. 
We can then use the Fisher-Hartwig conjecture [33] to rigorousljf|] calculate the large 
A'^ behavior of Sa- However, going through all the technical complications of the 
Fisher-Hartwig conjecture is not needed, because we can exploit the results obtained 
for the very similar matrices of lattice free fermions on the infinite line. In fact, for 
lattice models in the thermodynamic limit at fixed filling v = A^/Liat, it has been 



found that the entanglement entropies of £ut consecutive sites is given by Eq. (21) 



§ Despite its name, the Fisher-Hartwig conjecture has been rigorously proven for the case at hand 
(see e.g. Ref. l34l '). 
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2.5 1 ^^ 



(Z3 




3 4 

In [2N sin(7U X)] 



Figure 1. Sa(N) for a = 1,2,5, oo (from top to bottom) as function 
of ln[2Af sin(7rX)] with X = i/L calculated numerically for several values of 
N = 10, 20, 40, 80, 160. The straight lines are the asymptotic predictions in Eq. 
| |31[ l. The convergence is clear, although for a > 1 non-monotonic corrections to 
the asymptotic behavior are present. 



where 15^4 (A) is a standard tut x ^lat determinant with the correlation matrix C given 
by [3011^ 

sin7rz^(?i — m) 



^lat 



■K[n — m) 



n.m ■■ 



1,...,^] 



lat 



(30) 



It is evident that this matrix C^^ is the same as A in Eq. (29) identifying i/ with 
i/L. However, this is only a mathematical coincidence and it will most probably not 
be true for interacting systems. Indeed, in Eq. (29) we have a finite system and the 
indices are related to occupation modes of the N particles in the full system. On the 



other hand, in Eq. (30) we have an infinite lattice with filling ly and the indices refer 
to the lattice sites of the subsystem. Having established this equivalence between the 
two matrices, we can use the exact calculations in Ref. [3S] (see also [5T]), replace v 
with i/L, and obtain the asymptotic behavior of the desired entanglement entropies 
as 



On 



1 



1 



L 



In ( 2N sin tt-] +Eo,+ o{N"). 



(31) 



Eq. (|3lj) agrees with the CFT prediction for finite systems with periodic BC in Eq. 
(Isl) and represents an explicit analytic confirmation of this CFT result. It coincides 
with the scaling prediction in Eq. ( 13 ) from the lattice model, but here it has been 
derived from first principles. 



Notice that we cannot recover the infinite volume limit from Eq. ( 29 ) because 
this limit must be taken at fixed ratio N/L. If we naively take L — >■ oo in Eq. (29) 
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8 0.01 




0.001 r 



100 



Figure 2. Leading asymptotic correction to the Renyi entanglement entropies of 
half system i = L/2 for N up to 100 and a = 2, 5. Straight lines correspond to 
the prediction | |33| l that agrees with numerical data. Note that also subleading 
corrections oscillate and can be described by Eq. ||34|. 



we get a meaningless result, reflecting the non-commutation of the limits. Oppositely, 
after computing the determinant as in Eq. (31), the infinite volume limit exists at 
finite density N /L. 

Fig. [T] shows a comparison with exact finite- iV calculations, for a = 1,2, 5, oo. It 
is evident that (especially for large a) the data are affected by finite N corrections 
that are exactly calculated in the next subsection. 



3.2. Corrections to the asymptotic behavior and universal FSS in finite chains 

The above correspondence, between the determinants giving the entanglement 
entropies of the continuos system and the ones for spin chains, permits a quantitative 



description of the scaling correction to the leading behavior in Eq. (31 ) by exploiting 

[201 El]. 



the results based on the generalized Fisher-Hartwig conjecture [36j in Refs. 

We introduce the differences between the entanglement at finite N and the asymptotic 



values S^^y{N) in Eq. ^ as 

d^{N) = Sa^iN) ~ S^'^N) . (32) 

We can again use the spin-chain results of Refs. [20l [21], where a quantity 
analogous to da{N) was calculated at the leading order. Using these results and 
replacing ly with i/L (in [5H1 [SUl [5T] kp = tti^ was used, but here we prefer to use i/ to 
avoid confusion with the Fermi momentum in the continuous system) , we obtain that 
the leading correction term is given by 



do.{N) 



2cos{2nm/L) 
l-a 



(2A^sin7r^/L) 



-2/q 



^ V2 ^ 2a) 



0{N- 



i[4/q,2] 



(33) 



A check of the correctness of this expression is reported in Fig. [2| where we report the 
absolute value of da{N) for the half system entanglement (i.e. £/L — 1/2) for N up 
to 100. The power law behavior is evident and the straight lines are given by Eq. (33) 
without any adjustable parameter. These corrections of the form A^"^/" correspond to 
the ^~2/a corrections found within conformal field theory [3^, that have already been 
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generalized to other situations, such as massive field theories |38], confined systems 
[55] . disordered models |3D] and have been carefully checked numerically in many 
different models [201 US [12] . 

Subleading corrections to the scaling are visible in Fig. [2] and for large values of 
a they have a sizable effect. These can be exactly calculated adapting the results of 
Ref. [21] (based on generalized Fisher-Hartwig [36] and random matrix theory [42] ) 
and by replacing v with t/L. The full result for da{(-) up to order N~^ can be cast in 
the form 



2 "~ _ 2 



2p(2g-l) 



{Q.r 



p,q=l 



1 a + l 
Z|^285a3 



cos{2pttN£/L) Aq sm{2pTrNe/L) 



L 



N 



[Bp^qC^'P''^'^/^ + h.c. 



{ibiia^ - 7) + (49 - q2) sin2(7r£/L)) +o{l^' 



where 



Ln = 2A^sin(7r£/L), 



A, = 



Qq = 



Bn 



/ 
1 + 3 
\ 


2q-l 
a 


\n\ + '-^)] 


[n\-\^)\ 


2q-l 


r.^ + 7 



cos{n£/L) 



(34) 

(35) 
(36) 

(37) 



6a 

_ P 
4 



'(2'^~l)0sin^(-^A)-15(^^^^ + l) 



,(2g 






a^ 



cos{tt£/L) 



(38) 



As for the spin chain, while Eq. (34) provides an infinite number of contributions, 



for a given fixed value of a only a finite number of them will be smaller than the 
leading neglected term, which is always of order 0{N~^). To be specific, in the cases 
a = 2,3,10 Eq. (34) gives the leading 4,8,46 terms in the asymptotic expansion of 
da{N) and hence the leading 6, 10,48 terms in the expansion of Sa(i)- We do not 
report here all the terms which contribute for specific values of a that can be obtained 
from a simple adaptation of the results above or from Ref. [2T] . 

A remarkable exact result that we obtain from the previous analysis is the 
universal finite size scaling (FSS) form for finite XX spin chains with periodic BC. 
Indeed, as explained in section [Ll] the above result represents the continuum limit of 
the XX spin chain in a finite volume. For these spin chains, a universal FSS form has 
been observed in Ref. [20]. The quantity considered for the spin chain in Ref. [20] is 



Fc^iX) 



^— (^sin(7rX)2sin7ri^) 



-2/a 



r^ + s 



with X = Aat/ilat . (39) 



_^ V 2 2a > . 

In Ref. [20j chains with an odd number of spins iiat have been considered. Fig. p^ 
shows these results for a = 2 and magnetic field h — Q. The figure shows a perfect 
data collapse and somehow the correctness of the FSS ansatz. Already in Ref. [2D]j it 
was observed that the scaling function is perfectly described by F2{X) = ±cos(7rX), 
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Figure 3. Universal function for the corr ection to the scahng of the XX spin chain 
(lattice free fermion) _Fc«(^) in Eq. (J391 at zero magnetic field. Left: Numerical 
results for spin chains with o = 2, an odd number of spins 17 < L < 4623, and all 
possible values of ^i^t < -Liat/2. Data show perfect data collapse to the function 
-F2(^) = i costtX. Right: numerical results for Li^t = 101 against our prediction 
F^iiX) = cos(27r!/£iat), showing graphically that the curve in the right panel is 
just the envelope obtained by sampling the oscillating function at integer values 
of £i^i (see the text for details). 



as evident result from the figure (where it is impossible to distinguish the data from 
the conjecture). However, this is in apparent contradiction with our result, suggesting 



that a FSS form does not exist and the quantity Fa{X) in Eq. (391 should instead be 
just a function of iut and in particular 



F„( 



-latj 



cos(27ri/£iat) , 



(40) 



where we changed the variable of Fa from X to £iat . To elucidate what is happening in 
the right panel of Fig. [3] we report the numerical data at Liat = 101 in zero magnetic 
field against our new prediction Eq. ( 40 ) . 

Being Ljat odd, the ground state is not exactly at half-filling but at i' — 
{Liat — l)/2-^iat- In Fig- [3] (right panel) we can observe that the prediction (40), i.e. 



the strongly oscillating curve, agrees with the numerical data, apart some subleading 
corrections to the scaling. The ±cos(7rX) form (also shown in the right panel) is 
nothing but the sampling of the curve at the integer values of £iat . Indeed from basic 
trigonometry we have 



cos(27rj/£ 



lat ) 



(' 



-^lat 



1 



cos TT^lat cos TT 



Clat 



t;iat 



(TT^lat 



■ sin TT^iat sin tt 



Lla.t'' 

Aat 



(41) 



^lat 



(-l)^"''cos7rX, 



ilat 

where we used that ^lat is an integer. The final expression is exactly the 
phenomenological result conjectured in Ref. [2D] that we then prove and 
generalize here. Indeed for even chains, always in zero field, Ref. [20j proposed 
phenomenologically Fa{X) = ±1, which corresponds to cos(27ri/^iat) i.e., for v = 1/2, 
to (—1)^'"'. We checked Eq. (40) against other spin-chain results at different filling. 



always finding agreement. Leading corrections to the scaling having a structure similar 
to Eq. (40) have been conjectured in Ref. |43j for spin chains with open boundary 



conditions. 
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Figure 4. Universal function Foo {X) for the correction to the scaling of the single 
copy entanglement in the XX spin chain at zero-magnetic field. Left: Numerical 
results for spin chains with an odd number of spins 17 < L < 4623 and all possible 
values of £iat < Z/iat/2. Data shows perfect data collapse to the function -Foo(^) 
in Eq. Esl shown as a continuous line, indistinguishable from data poi nt. Right: 
numerical results for Li^t = 37 against our prediction FaaiX) in Eq. (J44|. It is 
graphically evident that the curve on the left is the envelope obtained by sampling 
the oscillating functions at integer values of ^i^t- 



3.2.1. The single copy entanglement: a —>■ oo. For a — >■ oo, the Renyi entanglement 
entropy gives the logarithm of the maximum eigenvalue of pA also known as single copy 
entanglement E|. It is not possible to obtain the result at a ^- cx) from the general 
form in Eq. ( |34[ ) because all the corrections of the form N"^^" resum. Again instead of 
re-doing all the calculation to resum these terms, we exploit the correspondence with 
the infinite spin chain, and simply obtain the final result from Ref. |21j substituting 
1/ with £/L. After straightforward algebra we obtain 

1 \u,{~e^'-'^/^)+Lh{^ 



dUN) 



-i27TiN/L\ 



(42) 



UnibNsiniTrl/N)] 

with b — exp(— ^(1/2)) w 7.12429. We have checked these results against exact 
numerical computation that we do not report here. 

It is interesting also in this case to explore the consequences of this result for 
finite spin chains, on the same lines as above for the finite a results. In Ref. |20| . 
for chains with an odd number of spins, it has been shown that the data for several 
choices of iiat and ^lat collapse on a single master curve if plotted as 



Foo{X) 



tlatj-t^lat) 



In ( 2b^:^^ sin(7rX) sin{Tn/) 



withX = -^lat/ilat ■ 



(43) 



Numerical data showing this collapse (analogous to those in Ref. i20j) are reported 
in Fig. [4] (left) for zero magnetic field and odd Liat- Oppositely to the case for finite 
a, the shape of this curve was too complicated to be guessed in Ref. [2^. As before, 
assuming universality in the FSS towards the continuum limit, we predict 



Foo(^lat) = ^ [Li2(-e' 



i2iny£i^t 



)+Li2(-e 



— i27rl/€lat^ 



(44) 



For a small chain with iiat = 37, the above curve is reported in the right panel of Fig. 
|4]and it shows high frequency oscillations, but perfectly coincides with the exact lattice 



The entanglement entropy of ID systems in continuous and homogenous space 



13 



0.01 


^ 




g 


= 









0.0001 




• 
■ 
* 


X=l/2 
X=l/3 
X=l/6 


NN^ 




\ 
















10 










-0.1 



-0.2, 



' 1 ' 1 ' 1 


^^ 


\I^^^y^ 


/^ 


\f ^T^ ^r 




ID jr^ ^r 




'_^^0^ X 




(i y 


• v=l/2 


L >^ 


■ v=l/3 
• v=l/4 - 


-.'i 1 1 1 1 1 


' v=l/5 


1 1 1 



0.1 



N 



0.2 0.3 
X 



0.4 



Figure 5. Left: Leading correction to the scaling for the von Neumann 
entanglement entropy for the free-fermion gas for several length of the subinterval 
£ = XL as function of A'^ up to 60. The straight lines are the asymptotic prediction 
in Eq. (Um. Right: Universal function Fi(X) in Eq. | |48[ l for the correction to 
the scaling of the entanglement entropy Si in the XX spin chain at zero-magnetic 
field as function of X = ^lat/^iat- Data are reported for two values of Li^t 
of order of 1000 and 300 which are compatible with the corresponding filling 
u = 1/2, 1/3, 1/4, 1/5. The continuous lines (indistinguishable from data points 
for not too small X) are the conjecture in Eq. |49|. 



calculations (apart small subleading corrections). As for finite a, the smooth result 
obtained for chains of different length (reported in the left panel) is a consequences of 
the sampling at integer fiat- Using the property of the Li2(y) function, the envelope 
in the left panel of Fig. [4] is 



Fao{X) 



[Li2(- 



ilTX\ 



Li2(-e 



-i-ITX\ 



(45) 



that is also shown in both panels, but in the left one is indistinguishable from data 
points. 



3.2.2. 



The von Neumann entanglement entropy: a = 1. For a = 1, the leading 
correction in Eq. ( 33 1 is vanishing. The actual calculation requires a complicated 
mapping with random matrix theory, but the final result, correct at 0{N^^), can be 
read from Eq. (34) specialized to a = 1: 

1 



diiN) = - 



127V2 V5 



+ (cot7r£/L)^ = - 



1 



12iV2 



1 



(46) 



Fig. ^ (left panel) reports the exact numerical computation for the entanglement 
entropy for various values of i/N, showing perfect agreement with the asymptotic 
formula. 

We consider again the XX spin chain for which the corrections to the scaling for 
the von Neumann entropy have not been considered quantitatively in finite systems. 
For infinite systems Ref. [21] reports the exact result 



dl(^lat)=- 



1 



12C 



(cot TTI^) 



(47) 
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It should be pointed out that this term is an "analytical correction" to the scaling, i.e. 
it is not due to the insertion of an irrelevant operator. For this reason, its finite-size 
scaling cannot be obtained simply by replacing the distance £ia± with the chord length, 
as done for da at finite N. 

For finite systems, we expect the FSS form 

dl(^lat, ilat) = ^ -2^l(^lat/^lat) , (48) 

(^Sin(7r£iat/£lat)) 

where Fi{X) is an unknown function with -F'i(O) fixed by Eq. ( [47| ). However, by 
looking at its continuum limit, it is reasonable to propose the FSS ansatz Fi{X) = 
A + Bsin'^ ttX. The constant A can be fixed by requiring that -F'i(O) is given by Eq. 
( [47| . The constant B can be fixed with the numerical data (e.g. by the scaling at 
X — 1/2). After a careful analysis, we conjecture the FSS scaling function 

Fi{X) = ( - +cot2 7ri^) H l + -cot2 7ri/ sin^iirX). (49) 

12 \ 5 / 18 5 



By construction, this form reproduces Eq. (47 1 for X = 0. In the continuum limit, i.e 



i^ — >■ 0, it reproduces Eq. (46), as an highly non-trivial check. A stringent test of its 
correctness is provided by numerical data. These are reported for few different values 
of I' and for large values of Liat in Fig. [5] right panel, showing a perfect agreement. 

4. Excited states in periodic chains 

We now turn our attention to excited states that can be easily treated in the formalism 



we introduced. Indeed, the only change compared to the ground state is in Eq. ( 16 1 
where we have to sum over the occupied one-particle levels. It is convenient to have a 
simple graphical representation of the many-body states. This can be easily done by 
representing each single particle state with a circle and filling in black the occupied 
ones and leaving empty the others. For example, the ground state is 

•••00«---«««---«00--- (50) 

"^ V ' 

N 

where the underlined circle represents the zero momentum mode. When working at 
fixed number of particles TV, excited states are obtained from the ground states just 
by moving black circles to empty white ones. 

The entanglement of excited states has been already considered few times in the 
literature, but only in the context of discrete lattice models. In [H] it was shown that 
the negativity (which is a different measure of entanglement, related to some Renyi 
entropies fV) shows a universal scaling in critical spin systems. In Ref. |45| . on the 
bases of Toeplitz matrix arguments for the XX spin chain and by exact calculations 
for the anisotropic Heisenberg one, it has been shown that only a small subclass of 
excited states can exhibit a universal logarithmic divergence with the subsystem size 
i, while most of the states strongly violates the area law and their entanglement 
entropies increase linearly with £, with non-universal prefactor. The states providing 
universal scaling are those where there is a finite (and possibly small) number of sets 
of one-particle states occupied sequentially, as e.g. 

•••00«««000«««««00--- (51) 
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and the locations of these blocks of states is not essential. This set of states includes all 
low-lying excited states. In Ref. 53 it has been shown that the entanglement entropies 
of low-lying excited states display a universal finite size scaling that is different from 
the one in the ground state of Eq. (Is]). These can be calculated by means of CFT, 
because low-lying states in CFT language are described by the action of a scaling 
operator on the ground state. The states that are obtained by applying a primary 
operator to the ground states are of particular importance. In this case, the Renyi 
entropies for integer a have been related to the correlation functions of these operators 
in a a-sheeted Riemann surface [35] ■ We remand the interested reader to the original 
reference ^46j and we limit to quote the main result 

TrpS = [Ttp"a](gs)F[''\i/L) = c„ (^- sin(7r^/L) j F^^\i/L) . (52) 

where [TrpJ|](G5) is the ground-state value. F~^ [X] is the universal scaling function 
depending on the operator T whose action on the ground state gives the desired excited 
state. In particular F^ (0) = !> i-e., in the thermodynamic limit, all these low-lying 
states have entropies degenerate with the ground-state, in agreement with Ref. [i5] . 
Two sets of primary operators can be easily treated for a free boson compactified 
on a circle, which describes the thermodynamic limit of the free-fermion gas we are 
considering. First, the vertex operators V{x) for which Ref. [15] reports Fy {X) = 1 
(i.e., the entanglement entropies are the same as in the ground-state). In the free- 
fermion gas, this corresponds to the excited states obtained by shifting the ground- 



state (50 1 in momentum space, i.e. replacing all ki with ki + M with M arbitrary 



integer number. The matrix A is always given by Eq. (28 1, that depends only of the 



differences between the various momenta, and so it is exactly equal to the ground-state 
one, confirming the prediction Fy (X) — 1. 

The other operator considered in Ref. [46] is T = id(j), which has been found to 
have a non-trivial scaling function given by 

F!f\x)^{-l)°'i-sm{TrX)] detH, (53) 

for integer a. Here H is a 2a x 2a matrix with elements 

ifJT^fc \ tt{2j ~2 + x)/a ii j < a 




hk = < e*^^ - e^^fc and z, = <^ ' . (54) 

' ^ if^-^fc ' \7T{2j-2~x)/a iij>a ^ ' 

For a = 2, this reduces to the simple expression 

^(.)(^)^(7 + cos(2^X))^ (55) 

^ ^ 64 ^ ' 

but for any other a > 2 the explicit formulas are too cumbersome to be reported in 
their full glory. It must be mentioned that the analytic continuation of F^ {X) is not 
yet known, and so also the von Neumann entanglement entropy of this excited state 
is still unknown. However for small X it has been found 

4") {X) = l- ^ (a - a-i) + 0(X3) , (56) 

whose analytic continuation is obvious. 
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Figure 6. Scaling functions for the entanglement Renyi entropies in the lowest 
energy particle-hole excited state of a periodic system in Eq. l |57| . Left top 
panels report the scaling function F^"' (X) for a = 2, 3 as function of X, showing 
good agreement with the CFT prediction ) |53| shown as continuous lines. The 
lower left panel reports the scaling function for the von Neumann entropy that 
for small X agrees with the general expansion l |56[ l, reported as a dashed line. On 
the right the universal corrections to the scaling are reported for the half-system 
entanglement showing the behavior N~^'°' for the leading corrections. For a = 1, 
the corrections are monotonic (last panel) and effectively described by Eq. (|62[|. 



The excited state generated by the action of id<j) on the ground state is the 
particle-hole excitation obtained by moving one particle from the highest occupied 
level to the first available one, i.e., graphically 

(57) 



•00««' 



'•••' 



■•0«00' 



Af-1 



The corresponding N x N matrix A has then the first n — 1 rows and columns identical 



to the ground state Eq. (29), but the last different, given by 

sin7r(Af + 1 -■m)l/L 



7r(7V + l -to) 



1, 



,iV-l 



(58) 



l/L 



N 



and A„iAr — Aj^m- Despite only one row and one column differ from the ground-state, 
the matrix A ceases to be a Tocplitz matrix and (to the best of our knowledge) no 
analytic treatment is possible anymore. 
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We check the prediction in Eq. ( 53 ) numerically. In Fig. ml we report the 
numerical calculated scaling function 

F(")(X) ee exp[(l - a){SaiN) - (1 + l/a)/6lii{2N simrX) - ^„)] , (59) 

for several values of A^ as a function of X for a = 2,3. It is evident that in the large N 



limit the CFT prediction ( 53 ) is approached with small oscillating corrections to the 
scaling which are more pronounced for small X. In order to shed some light on the 
analytic continuation at a — >■ 1, we also report (always in Fig. pi) the scaling function 
for the von Neumann entropy 

F^^ (X) = Si{N) - 1/3 ln(2 AT sin t:X)~Ei. (60) 

As a difference with F^") with a > 2, the corrections to the scaling are much 
smaller, as for the entanglement in the ground-state. Unfortunately, as already 



stated, the analytic continuation to a —> 1 of Eq. (53 1 for arbitrary X is not yet 



known and so the data in the figure cannot be contrasted to an exact prediction. 



However, such an analytic continuation is known for small X: from Eq. (56) we have 
F^^(X) = 2(7rX)2/3 + 0{X^). This prediction is reported on top of the numerical 
data and they agree perfectly up to X ^ 0.1. 

Having established the leading asymptotic behavior, we move our interest to the 
leading corrections to the scaling. We find numerically that the corrections have the 
same exponents as in the ground-state, i.e. they decay with a dependent power-law 
^-2/a jj-^ order to show this, we report in Fig. ^ (left panels) the quantity 

Da{N) = [S^{N) ~ S^'^N)]N^^" , (61) 

for the half-system entanglement. The data clearly show the behavior for the 
corrections of the form (— l)^iV-2/a j^j, a > 1. Wc found numerically that 
-D2 -- (-1)^0.19039 ... and D3 - (-1)^0.225 . . .. These non-universal amplitude are 
different from the ones found for the ground state and we have been not able neither 
to calculate nor to guess their a dependence. We checked that for general £/L, the 
corrections are of the form cos{2itN£/L)N~'^^°' as for the ground state. Furthermore 
subleading corrections seem to have the same power structure as in the ground-state 



(cf. Eq. ([M])). 

The von Neumann entropy at a = 1 requires a separate analysis. Indeed, as the 
last panel of Fig. [6]shows, the correction to the scaling are monotonic. However in this 
case we do not know exactly the constant term in the leading behavior. An accurate 
numerical analysis for the half-system entanglement is consistent with the behavior 

5i(iV) = ^log27V + yi + ^, (62) 

with 2/1 — 0.540726 . . . and 2/2 = 0.35 These numerical data have been obtained 

by fitting data for N > 100, 150, 200 and keeping under control the stability of the fit. 
Although these fitting parameters have been extracted from asymptotic large N, Fig. 
|6] shows that the fit describes very accurately the data down to A^ ~ 3. 

5. Systems with boundaries 

5.1. Hard-wall boundaries 

We now consider a gas of spinless fermions confined in the interval [0, L] by a hard- 
wall potential, i.e. the gas density vanishes outside the interval x ^ [0, L] and the 
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Figure 7. On the left we report the Renyi entanglement entropies 5a for the 
free-fermion gas confined by hard wall potential for o = 1,2,3 (from top to 
bottom) as function of ln{4A'sin7rX) with X = tjL. For the larger values of a, 
osci llating corrections to the scaling obscure the leading asymptotic behavior Eq. 
(J65l represented as straight lines. These oscillating corrections are quantitatively 
described by Eq. l |67[ l that is checked in the right panel for the same values of a 
for the half-system entanglement (£ = L/2). 



boundary condition is that the wave-function vanishes at the boundaries (Dirichlet 
BC). The one-particle wave functions are 



with energies E^ — 7r^fc^/2_L^ 



2 

— sm 
L 



1 ^ 
Trfc — 

L 



1,2,.. 



(63) 



5.1.1. An interval starting from the boundary. The elements of the overlap matrix 



(cf. Eq. ( 23 )) between two one-particle eigenstates n and m have a particularly simple 
form for an interval starting from the boundaries, i.e. A — [0,^]. In fact we have 



A„„i = / dz (j)*^{z)(t)m{z) = 

Jo 



(P\ = sin[7r(n - 
u*; — 



m)e/L] sin[7r(n -|- m)i/L] 



(64) 



■n(n — m) 7r(n + m) 

with n,m = 1, ...,N . 

As for the periodic case, the matrix A above is exactly the same as the correlation 
matrix C^''* of an infinite lattice with v replaced by (^/N . This has been considered 
in Ref. [35] where, using a recent generalization of the Fisher-Hartwig conjecture to 
Toeplitz+Henkel matrices [17] , the asymptotic behavior of the entanglement entropies 
for the lattice model have been calculated exactly. Exploiting the equivalence between 
the two problems (i.e. replacing v with £/iV in Ref. [13]) we easily obtain for the 
asymptotic behavior of the entanglement entropies 



bo. 



hi^ 



a.) 



1) sinTrX] 



Ea 

~2 



(65) 
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where Ea is defined in Eq. (11). Notice that this result agrees with the general CFT 
prediction in Eq. Q with Ing = that is a well-known result for open boundary 
conditions HI}]. 

A comparison of the finite- TV results with Eq. ( 65 ) is shown in Fig. [t] It is evident 
that for any a there are corrections to the scaling oscillating with N . These are of the 
order 0{N~^/°') and can be deduced exactly from the analogy with the lattice model 
solved in Ref. [43]. Defining 

d^{N) ^ S^{N) - S^^^N) , (66) 

we have from Ref. [13] and replacing v with t/L 

1-" r(2-2^) 

Fig. [t] (right panel) show these corrections for half-system entanglement entropy 
for a = 1,2,3. Further corrections of the form N~P'°' with p integer can be 
straightforwardly deduced from the analysis in Ref. [43 . We mention that, as for 
the periodic case, these leading N~^/°^ corrections correspond to the ones of the form 
0(^-i/") found within CFT [371113]. 

5.1.2. Generic interval. In the case of a subsystem consisting of generic interval 
A — [xi,X2], the entanglement entropies require a different analysis. The general 
formula for the matrix A is slightly more complicated: 

A„m = / dz(t)*^{z)(j)m{z) =^nm{x2)-^nrn{xi), (68) 

J X\ 



where the matrix B is defined in Eq. (64). The entanglement entropies of the interval 



[a;i,a;2] can be computed by inserting its eigenvalues in Eq. (26). This allows us to 
easily compute Sa{N) up to large values of TV, and compare its behavior with the 
asymptotic CFT prediction 



K-^){ 



1 



6V a J ■ " ■ 2 



sin^[7r(y2 - yi)/2] sin(7ryi) sin(7ry2) 



S„(iV) = - l + -Uln47V + S„ + -ln ' '"^ " ' ' '-^^ '-^^ ^ (69) 



sin^[7r(2/2+2/i)/2] 



where yi = Xi/L. The proof of this equation is a straightforward CFT exercise that 
we report in [Appendix A[ 

In particular, we considered a block of size Ly centered at the middle of the 
system, i.e., 2/2 = 1/2 -|- y/2 and j/i = 1/2 — y/2. The data of the entanglement 



entropy approach the asymptotic behavior predicted by Eq. (69 1, i.e., 

Sa,iN) = l(l + -) [In 2A^ + ^, + In sin(7ry)] , (70) 

6 V a/ 

as shown for a = 1 in Figs. Is] (left panel). It can be seen numerically that the 

corrections to the scaling are of the order 0{N~^^°'). It must be noted however that 

the 0{N~^'") convergence is nonuniform approaching the limits x —>■ and a; — > 1, 

where the entanglement entropy trivially vanishes for any N. Simpler results are 

obtained also considering an interval starting from the center, i.e. taking X2 — 1/2 ±y 

and Xi — 1/2, for which 



SaiN) = ^(1 + ^) |ln4TV + £;„ + ^ In [cos(7r2/) tan(^2;/2)^ 



,(71) 



see Fig. [s] (right panel) for a = 1. We also checked the correctness of Eq. ( |69| for 
other values of a and for different choices of xi and X2, but the resulting figures are 
not very illuminating and we do not report them. 
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Figure 8. (Color online) The von Neumann entanglement entropies Si of 
intervals [a;i,a:2]- In particular we consider the cases x\ = L{l/2 — y/2), 
X2 = L{l/2 + y/2) (left) and xi = L/2, X2 = L{l/2 + y/2) (right). We plot 
Si (Af) — (In N)/6 against y. In both figures the black full line shows the asymptotic 
behavior given by Eq. l|69|. 



5.2. Neumann boundary conditions 

Another interesting situation arise when imposing Neumann boundary conditions on 
the fermionic wave-function, i.e. imposing that the derivative of the wave function 
vanishes at the two boundaries at and L. In this case, the normahzed one-particle 
wave functions are 



Mx) = 



2-S, 
L 



fe.O 



cos 



nkx 



fc^O, 1, 



(72) 



with the same energy as for Dirichlct BC. As an important difference compared to 
Dirichlet BC, also the zero-mode with fc = does contribute. 

The A matrix is readily calculated. It is an N-hy-N matrix with entries that are 
more easily written if we count rows and columns with n, m starting from and up to 
TV — 1 as for the modes above. For an interval of length t starting from the boundary, 
straightforward calculations lead to 



A, 



sin[7r(n — m)z\ sin[7r(n -|- m)z\ 



7r(n — m) 



ii{n + m) 



if n, TO = 1, 



.A^-1, 



72^ 



if n = and m ^ 0, 
ii m,n = 0, 



(73) 



and Aom = A„jo- Note the plus sign between the two terms for n,m ^ Q and the 
zero-mode contribution, as a difference compared to Dirichlet BC. 

Because of the presence of the zero row and column, A is not of the form 
Toeplitz+Hankel as it is for Dirichlet BC. Thus the recent generalizations of Fisher- 
Hartwig conjecture in Ref. [47] cannot be used. We then determine numerically the 



matrix A for various N and, through Eq. ( 26 ) , we compute the Renyi entanglement 
entropy shown in Fig. ^ (left panel). The analysis of their large- A'^ behavior gives 



So.{N) 



(74) 



l^{l + \) ln[2(2Ar - 1) sin(7r^/i)] + ^ + • • • , 

shown as continuous lines in the figure. This form is consistent with the general CFT 
expectation in Eq. (H) with g = 1. In order to avoid confusion with CFT literature. 
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Figure 9. On the left we report the Renyi entanglement entropies Sa for the 
frec-fermion gas with Neumann BC for a = 1,2 (from top to bottom) as function 
of ln(4Af siuTrX) with X = (./ L. Oscillating corrections to the scaling are evident 
on top of the leading asymptotic behavior Eq. ( |74[ l represented as straight lines. 
These oscillating corrections are quantitatively described by Eq. | |75| l that is 
checked in the right panel, where we report the quantity r^ in Eq. l |76| for 
ct = 1,2,3 against the theoretical prediction from Eq. l |75[ l for the half-system 
entanglement {£ = L/2). 



we stress that, in this paper, we are considering Neumann and Dirichlet BC on the 
ferniion degrees of freedom. These do not correspond to Dirichlet and Neumann BC on 
the bosonic field obtained from the bosonization of the fermionic theory that instead 
are well known to have different g function (see e.g. [48 ). It is known that they both 
correspond to Neumann conditions of the bosonized field and so it should be not a 
surprise that the asymptotic behavior up to 0{N'^) is the same as Eq. (65 1. Notice 



that we have included a 0{\/N) term in the leading behavior of the logarithm (i.e. 
the —1 in (27V— 1)) that has the effect to cancel the leading non-oscillating correction 
to the scaling. This was present also for Dirichlet BC, but it has opposite sign. While 
before this was motivated by the mapping to the lattice model (cf. Ref. jlS]), here we 
introduced it on a phenomenological basis in order to describe the data (see below) 
and we do not have any mathematical explanation for it. 

We now consider the corrections in N to the leading behavior, that are again 
consistent with the general scaling from CFT 0(7V~^/"). On the basis of the numerics, 
we guess exactly the first correction to the scaling, and we can write the entanglement 
entropies as 

ScW = 1 f 1 + i) \n[2{2N - 1) sin(7r^/L)] + ^ + 



12 V"^ ■ a. 

2sin[7r(2Af ■ 



+ 



l)i/L] 



[2(2N ~ I)sin7r£/L] 



-1/' 



! V2 ^ 2a I 



2a I 



(75) 



In Fig. [9] (right panel) we show the evidences for this scaling for the half-system 
entanglement (i = i/2). We report the quantity 

_ \S^{N) - \ ln(2(2iV - 1)) - EJ2\ 



[2(2iV- 1) 



-1/q 



(76) 
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that for a = 1,2,3 approaches for large N the value predicted from the ansatz (75) 
fa = Tzr^ ril-AL\ - Fo'^ a = 1 the leading corrections are of order 1/A^. The figure 
shows that the absolute value for even and odd N coincide (but they have opposite 
signs). This confirms that choosing to parametrize the leading term with (2A'^ — 1) in 
Eq. (|74|) cancels the 1/A^ non-oscillating corrections completely. Some of the factors 



{2N— 1) in Eq. ( 75 ) are subleading and are not explicitly tested by the numerical data 



presented. They have been introduced from the analogy with Dirichlet BC results. 



Finally, we stress that we do not have any mathematical basis to justify Eq. (75): 
while its general structure can be inferred by CFT [37] because fermionic Neuman BC 
are in the same universality class as Dirichlet BC, the non-universal amplitude of the 
correction has been guessed and conjectured by exploiting the analogy with Dirichlet 
BC and tested agains numerical data. 

We also considered numerically other situations, such as other bipartitions of the 
systems. However, none of these results present particularly relevant or unexpected 
features to be mentioned here. 

6. Conclusions 

In this manuscript we report the details about the computation of the entanglement 
entropies of continuous systems (gases) which have been anticipated in the short 
communication [TH]. The most important ingredient to write down the entanglement 
entropies in terms of finite determinants is the use of the reduced overlap matrix in Eq. 



( 23 ) . The calculation of the entanglement entropies is then mapped to the solution of 
an eigenvalue problem of an A^ x iV matrix, with N being the number of particles of 
the gas. 

For the ground state of a periodic system we obtain the leading behavior in the 
form 

Sc. = Ul + ^yn (2N sin TT^\ +E^+ o{N°), (77) 

while for a gas with Dirichlet or Neumann boundary conditions we find 

^. = ^ (1 + ^) In (^msiuT:^^ +^+ o(iV°), (78) 

both in agreement with CFT and scaling expectations, but they have been found 
here from first principles. We also derive the corresponding leading behavior for some 
classes of excited states. 

Furthermore, adapting to the problem at hand the results in Refs. [2I1I27]! we 
calculate also subleading corrections. The universality of these formulas allowed us to 



infer the finite-size scaling forms for spin chains which are reported in Eq. ( 40 ) for 



1 < a < 00, in Eq. (42) for a = 00, and in Eq. (49) for a = 1. The determination of 
these exact formulas were left as open problems from previous investigations. 

Some other applications (such as to systems with defects, star graphs, and to 
gases confined by an external potential both in and out of equilibrium) of this novel 
method have been already shortly presented in Ref. |T9], but they will be detailed 
elsewhere. Other generalizations, which we are currently investigating, concern the 
calculation of the entanglement for quadratic Hamiltonian which do not conserve the 
fermion number (such as the continuum limit of the XY model) , free gases in higher 
dimensions and different geometries. Finally, some non equilibrium situations such as 
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local quantum quenches (e.g., instantaneously turning on/off of a defect) can also be 
tackled within this framework. The asymptotic CFT results in several circumstances 
are known [ST], but analytic calculations for specific models are still missing. They may 
provide important insights in view of the recent proposals of using the full counting 
statistics after a quench as an experimental probe and a measure of entanglement [52] . 
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Appendix A. The CFT entanglement entropy of an arbitrary interval in a 
finite system with Dirichlet boundary conditions 



We provide here the CFT proof of Eq. ( 69 1 . In CFT the moments of the reduced 



density matrix can be written as correlation function of particular twist fields that 
transform as primary operators under a conformal transformation [6l El [49]. In 
particular, for the case of a finite interval between xi and X2 in a boundary theory, 
we have (for integer n) 

Trpl = ($„(xi)$_„(x2)), (A.l) 

where $„ and $_„ transform as primary fields with dimensions a;„ — (c/24)(n— ^/n). 
This two-point function in the finite strip with Re(w) G [0, L] can be obtained 
from its conformal mapping to the upper half plane (UHP) with Im(z) > 0. The 
mapping and its inverse are 

m=-lnz, z = e^™/-^. (A.2) 

Using the property that the twist fields for free fermions with open boundary 
conditions behaves like primary operators in a free bosonic theory [43j . the two-point 
function in the UHP can be read from Ref. [SD] 

($„(zi)<i>_„(z2))uHP = c„ { ^^^^"^ \ " , (A.3) 

\z12Z12Z11Z22 J 

with Zij — \zi — Zj\ and zj. = 'Zk and c„ an undetermined constant (we set the UV 
cut-off a to 1). We stress that in general (i.e. for a theory that is not free bosonic) a 
universal function of the harmonic ratio build with the four points Zi and Zi multliplies 
the above formula [3TJ [35] . Using then the conformal mapping (|A.2 ) , we have 



z,2zr2 = (e^™^/^ - e--^/^)(e-'™i/^ - g— ■"'^/^) = 4sin2 ^^""^ "^^^ , 

ZIj 

\zii\ =2sm-^, Izaal = 2sm-^, 

212^12 = (e*™i/-^ - e-''^'^2/-^)(e-*™i/-^ - e*™^/^) = 4sin2 ^(^1+^2) ^^^^^ 

ZIj 

Thus 



($„(r, r)$_„(0, r)),trip = ^'(^1)!""" ^'(22)!""" ($„(zi(u;))$_„(z2(u;))) 



UHP 



2 „:„2 7r{wi+W2) \ 2="=" 

2L (A.5) 



^^^ sm ^ ^^ ' sm 
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Using Eqs. (A.l) and (fTj), performing the analytic continuation from n to a, we 



have 



r, r -1 -2 7r(uji— u;2) • 7rw-[ 

2L 1 sm ^ ^r — sm ^^ sm 



7TW2 



6V aJl IT 2 -2 77(m.i+«»2) J «' V 



sm 2L 



where the constant Ea is given by Eq. (11), according to the result for the single 



interval in a periodic systems. Assuming now the scaling hypothesis when working 



with finite number of particles, Eq. (69 1 follows simply by replacing L/tt by 2N using 
the argument in Section [l.l[ 
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